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ABSTRACT 

Context. Deuterium fractionation is known to enhance the [DCO + ]/[HCO + ] abundance ratio over the D/H~10~ 5 
elemental ratio in the cold and dense gas typically found in pre-stellar cores. 

Aims. We report the first detection and mapping of very bright DCO + J=3-2 and J=2-l lines (3 and 4K respectively) 
towards the Horsehead photodissociation region (PDR) observed with the IRAM-30m telescope. The DCO + emission 
peaks close to the illuminated warm edge of the nebula (< 50" or ~ 0.1 pc away). 

Methods. Detailed nonlocal, non-LTE excitation and radiative transfer analyses have been used to determine the 
prevailing physical conditions and to estimate the DCO + and H 13 CO + abundances from their line intensities. 
Results. A large [DCO + ]/[HCO + ] abundance ratio (> 0.02) is inferred at the DCO + emission peak, a condensation 
shielded from the illuminating far-UV radiation field where the gas must be cold (10-20 K) and dense (> 2 x 10 5 cm" 3 ). 
DCO + is not detected in the warmer photodissociation front, implying a lower [DCO + ]/[HCO + ] ratio (< 10~ 3 ). 
Conclusions. According to our gas phase chemical predictions, such a high deuterium fractionation of HCO + can only 
be explained if the gas temperature is below 20 K, in good agreement with DCO + excitation calculations. 
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Molecules are enriched in d euterium over the elemental 
D/H abundance (1 - 2 x lQ- 5 . lLinskv et al.ll2006ft in many 
different astr ophysical environm ents. These include cold, 
dense cores dGuelin et al.l Il982ft. mid-planes of circum- 
stella r disks (jvan Dishoeck et all l2003t iGuilloteau et all 
120061) . hot molecular cores ([Hatched et al .11 19981 ). and even 
PDRs (jLeurini et alj 12006ft . Multiply deuterated species 
were first detected several years ago, e.g. D2 CO in warm 
gas (iTurnerl I1990D and NHD 2 in cold gas tiRoueff et al.l 
12000ft . iSolomon fc Woolfl l| 1973ft and IWatsonl (| 1974ft first 
proposed that deuterium fractionatio n is mostly caused 
by gas-phase ion-mol e cule re actions. ISmith et al.l {1982) 
and lRoberts fe Millarl (|2000al ) confirmed that the deuter- 
ation of H^~ at low temperatures (< 25 K) and of CH^ 
at higher temperatures (up to ~ 70 K) are important pre- 
cursor reactions in the subseque n t deuteration o f other 
speci es. iRoberts fc Millarl (|2000bft . I Walmslev et al.l (|2004ft 
and iFlower et al.l (|2006ft succeeded in reproducing the 
amount of several multiply deuterated molecules in cold 
gas by adding to pure gas-phase chemistry the accretion 
(freeze-out) of gas-phase molecules onto the surfaces of dust 
grains. Finally the observed deuterium fractionation in hot 
cores is thought to result from the liberation of deuterated 
molecules, trapped in ice mantles in the prestellar phase. 

Although it has been studied thoroughly for 30 years, 
deuterium chemistry is not yet fully understood. With 
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many chemical and physical processes competing for ef- 
ficient fractionation, models are easier to compare with 
observations for sources with well described physical con- 
ditions. In this letter, we report the detection of very 
bright DCO + lines in the Horsehead edge. In particu- 
lar, the mane of the Horsehead nebula is a PDR viewed 
nearly edge-on ( inclination < 5°) illuminated by th e 
Q9.5V star a O ri (lAbergel et al.ll2003HPhilipp et alll2006ft . 
lHabart et al.l (|2005ft showed that the PDR has a very 
steep gradient, rising to uh ~ 2 x 10 5 cm~ 3 in less than 
10" or 0.02 pc, at a roughly constant thermal pressure of 
~4x 10 6 Kcm~ 3 . The newly detected DCO + lines arise 
fro m a condensatio n adja cent to the PDR, first detected 
by iHilv-Blant et all ((2005). According to its submillimeter 
continuum emission, this core (B33-SMM1) is 0.13x0.31 pc 
long and has an average H2 densi ty of ~ 10 4 cm" 3 and 
a pea k density of ~ 6 x 10 5 cm -3 (jWard-Thompson et al.l 
12006ft . 



1. Observations and data reduction 

The DCO + J=3-2 line was observed during 2 hours of ex- 
cellent winter weather (~ 0.7 mm of water vapor) using the 
first polarization (i.e. nine of the eighteen available pixels) 
of the IRAM-30m/HERA single sideband multi-beam re- 
ceiver. We used the frequency-switched, on-the-fly observ- 
ing mode. We observed along and perpendicular to the di- 
rection of the exciting star in zigzags (i.e. ± the lambda 
and beta scanning direction). The multi-beam system was 
rotated by 9.6° with respect to the scanning direction. This 
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Table 1. Observation parameters. The projection center of all the data is a 2 ooo = 05 ?l 40 m 54.27 s , S2000 = -02°28'00". 



Molecule 


Transition 


Frequency 


Instrument 


# Pix. a 


TTi a 
^eff 




Resol. 


Resol. 


Int. Time ' 6 


Noise c 


Obs. date" 






GHz 










arcsec 


km s — 1 


hours 


K 




H 1;i CO+ 


J=3-2 


260.255339 


30m/HERA 


9 


0.90 


0.46 


13.5" 


0.20 


5.9/11.3 


0.06 


Mar. 2006 


H 13 CO+ 


J =1-0 


86.754288 


30m+PdBI 


2 


0.95 


0.78 


6.7" 


0.20 


2.6/4.5 


0.10 


Sep. 2006 


DCO+ 


J=3-2 


216.112582 


30m/HERA 


9 


0.90 


0.52 


11.4" 


0.11 


1.5/2.0 


0.10 


Mar. 2006 


DCO+ 


J=2-l 


144.077289 


30m/CD150 


2 


0.93 


0.69 


18.0" 


0.08 


5.9/8.7 


0.18 


Sep. 2006 


c 1B o 


J=2-l 


219.560319 


30m/HERA 


9 


0.91 


0.55 


11.2" 


0.11 




0.26 


May 2003 


Continuum at 1.2 mm 


30m/MAMBO 


117 






11.7" 











a Those columns apply to the 30m data but not to the PdBI data for the H 13 CO + J— 1-0 line. b Two values are given for the 



integration time: the on-source time and the telescope time. c Noise values estimated at the position of the DCO + peak. 
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Fig. 1. IRAM-30m integrated intensity maps. Maps have 
been rotated by 14° counter-clockwise around the projec- 
tion center, shown as the green cross at (Sx, Sy) — (20", 0"), 
to bring the exciting star direction in the horizontal direc- 
tion and the horizontal zero has been set at the PDR edge, 
delineated by the dashed blue vertical line. The synthesized 
beam is plotted in the bottom left corner. Values of contour 
levels are shown on each image lookup table. The emission 
of all lines is integrated between 10.1 and ll.lkms -1 . 



ensured Nyquist sampling between the rows except at the 
edges of the map. The DCO + J=2-l was observed dur- 
ing 11.3 hours using the C150 and D150 single-side band 
receivers of the IRAM-30m under ~ 8.5 mm of water va- 
por. We used the frequency-switched, on-the-fly observing 
mode over a 160" x 170" portion of the sky. Scanned lines 
and rows were separated by 8" ensuring Nyquist sampling. 
A detailed description of the C 18 J=2— 1 and 1.2 mm 
co ntinuum observat i ons a nd data reductions can be found 
in iHilv-Blant et all (|2005f) . We estimate the absolute posi- 
tion accuracy to be 3". 

We also use a small part of the H 13 CO + (J=l-0 and 
J=3-2) data, which were obtained with the IRAM PdBI 
and 30m telescopes. The whole data set will be compre- 
hensively described in a forthcoming paper studying the 
fractional ionization across the Horsehead edge (Hily-Blant 
et al. 2007, in prep). In short, the H 13 CO+ J=3-2 line was 
observed under averaged winter weather (~ 3.5 mm of wa- 
ter vapor) in rasters along the direction of the exciting star 
using the first polarization of the unrotated HERA. Each 




Fig. 

15". 



<5x (") 

2. Cut along the direction of the exciting star at Sy 



pointing of the rasters was observed in frequency-switched 
mode. This resulted in a 140" x 75" map, Nyquist sampled 
along the direction of the exciting star but slightly under- 
sampled in the orthogonal direction (i.e. rows separated by 
6" instead of 4.75"). The noise increases quickly at the map 
edges which were seen only by a fraction of the HERA pix- 
els. We finally used a frequency-switched, on-the-fly map of 
the H 13 CO+ J=l-0 line, obtained at the IRAM-30m using 
the A100 and B100 3mm receivers (~ 7 mm of water va- 
por) to produce the short-spacings needed to complement 
a 7-field mosaic acquired with the 6 PdBI antennae in the 
CD configuration (baseline lengths from 24 to 176 m). 

The d ata process ing was done with the GILDAsQ 
softwares (jPetvl l2005h . The IRAM-30m data were first 
calibrate d to the Tt scale us ing the chopper wheel 
method dPenzias &; Burrusl Il973f) . and finally converted 
to main beam temperatures (T m b) using the forward and 
main beam efficiencies (F c r & B c g) displayed in Table [TJ 
The resulting amplitude accuracy is ~ 10%. Frequency- 
switched spectra were folded using the standard shift-and- 



1 See http : //www . iram . f r/IRAMFR/GILDAS for more informa- 
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add method, after baseline subtraction. The resulting spec- 
tra were finally gridded through convolution by a Gaussian. 

2. Results and discussion 

Figure [Q presents the DCO+ J=2-l and J=3-2 and 
the C 18 J =2-1 integrated intensity maps, together with 
1.2 mm continuum emission. All maps are presented in a 
coordinate system adapted to the source geometry, as de- 
scribed in the figure caption. The DCO + emission is con- 
centrated in a narrow, arc-like structure, delineating the 
left edge of the dust continuum emission. A second max- 
imum is found at the extreme left of the map, associ- 
ated with a smaller dust continuum peak. Figure [2] shows 
the H 13 CO + and DCO + spectra in a cut along the direc- 
tion of the exciting star at Sy = 15" (horizontal dashed 
line of Figure [1]) . This cut intersects the DCO + emission 
peak which is close to the illuminated edge of the nebula 
(< 50" or ~ 0.1 pc). To our knowledge, this is the brightest 
DCO + emission (4K) detected in an in terstellar cloud close 
to a bright H 2 / PAH emitting region ([Habart et all l2005t 
iPetv et alj|2005h . The 15" spatial shift between the DCO+ 
and C 18 0/continuum emission peaks likely results from the 
steep thermal gradient. The region where the DCO + emis- 
sion is produced, is probably cooler than the region where 
the C 18 lines and 1.2 mm continuum emission peak, i.e. 
cooler than 30 K, the minimum temperature needed to ex- 
plain the intensity of the C 1 8 Q J=2-l lines in the cloud 
edge ([Goicoechea et al"]|2006f) . 

In order to constrain the [DCO + ]/[HCO + ] abun- 
danc^3 ratio from the observed line emission, we as- 
sumed that both species coexist within the same gas 
(implying the same physical conditions). This assump- 
tion is mainly justified by the spatial coincidence of 
the H 13 CO + and DCO + emission peaks (i.e. where 
T mh {DCO+(2 - 1)} /T mb {H 13 CO+(l - 0)} ~ 1). Besides, 
we used the H 13 CO + lines to determine the line-of-sight 
HCO + column density. Indeed, the direct determination of 
the HCO + column density from its rotational line emission 
is hampered by the large HCO + line opacities and their 
propen sity to suffer from self-absor ption and line scattering 
effects (|Cernicharo fc Guelinlll987f ). In addition, large criti- 
cal densities for HCO + (and its isotopologues) are expected 
even for the lowest- J transitions due to its high dipole mo- 
ment: ~ 4 D (Table [2]). Hence, thermalization will only 
occur at very high densities. For lower densities, n < n cr it, 
subthcrmal excitation dominates as J increases. Therefore, 
in order to accurately determine the mean physical condi- 
tions and the [DCO+]/[H 13 CO+] ratio at the DCO+ peak, 
we have used a nonlocal, non-LTE radiative transfer code 
including line trapping, collisional excitation and radia- 
tive excitation by cosmic bac kground photons (jGoicoecheal 
[200llGoicoechea et al.ll2006h . Collisional rates of H 13 CO+ 
and DCO + with H2 and He hav e been derived from the 
HCO+-H 2 rates of lFlowerl (|l999ft . 

Assuming a maximum extinction depth of 
Ay ~50 along the line -o f-sight where DCO + 
peaks ([Ward-Thompson et alj [2006), the observed 
DCO + J=2 1 and J=3-2 line intensities are well repro- 
duced (with line opacities ~1.5) only if the gas is cold 
(10-20 K) and dense (n(H 2 ) > 2 x 10 5 c m - 3 ). This high 
density is consistent with the one required to reproduce 



Table 2. Einstein coefficients, upper level energies and crit- 
ical densities for the range of temperatures considered in 
this work. 



Molecule 


Transition 








f'crit 






(s- 1 ) 




(K) 


(cm" 3 ) 


H i3 CO+ 


J=l-0 


3.9 x 10" 


-5 


4.2 


- 2 x 10 5 


H 13 CO + 


J=3-2 


1.3 x 10" 


-3 


25.0 


~ 3 x 10 6 


DCO+ 


3=2-1 


2.1 x 10" 


-4 


10.4 


~ 6 x 10 5 


DCO+ 


J=3-2 


7.7 x 10" 


-4 


20.7 


~ 2 x 10 6 



the CS J=5-4 excitation (|Goicoechea et al.l [2006) and 
with the value derived from dust s ubmm continuum emis- 
sion ([Ward-Thompson et al.l l2006h . The weakness of the 
H 13 CO+ J=3-2 line compared to the DCO+ J=3-2 line 
is caused in part by its larger Einstein coefficient (a factor 
~1.7 larger) and its higher energy level (see Tabled]). This 
implies that the H 13 CO+ J=3-2 line is more subthermally 
excited than the analogous DCO + line for the derived 
densities and temperatures. Note that we have not in- 
cluded collisions with electrons in this excitation analysis. 
In fact, the expected ionization fraction in su ch a cold and 
dense condensation is usually low, < 10 -7 ([Caselli et al.l 
I1999D . The derived DCO+ and H 13 CO+ column densities 
toward the DCO+ peak are ~ (0.5 - 1) x 10 13 cuT 2 
(i.e., [DCO+]~[H 13 CO+]~ (1 - 2) x lO" 10 ). Assuming a 
12 C/ 13 C=60 isotopic ratio ([Milam et al.ll2005l ). we finally 
find a [DCO+]/[HCO+] > 0.02 abundance ratio. 

In order to understand the observed deuterium frac- 
tionation in the dense gas close to the Horsehead PDR, 
we have modeled the steady state deuterium gas phase 
chemistry in a cloud with a proton density jih = "(H) + 
2n(H 2 ) = 4 x 10 5 cm" 3 illuminated by a FUV field 60 times 
the mean interstellar radiation field. We used the Meudon 
PDR codj^ L a photochemical mo d el of a unidimensiq nal 
PDR (see iLe Bourlot et al.l 119931; iLe Petit et ail l2006l for 
a detailed description) and its associated chemical reac- 
tion network. As this network only includes singly deuter- 
ated species, we added the D 2 and HP t s pecies and as- 
sociated reactions from iFlower et ail (|2()()(i ■. Nevertheless, 
these additional reactions do not affect much the pre- 
dicted DCO + abundances. Only H 2 , HD and D 2 form 
on grain surfaces because the used chemical network al- 
lows only H and D atoms to accrete onto dust grains. We 
chose the following gas phase abundances: D/H=1.6x 10 -5 , 
He/H=0.1, O/H=3xl0~ 4 , C/H= 1.4xlQ- 4 , N/H=8xlQ- 5 , 
N/H=8xl0- 5 , S/H=3.5xl0- 6 (| Goicoechea et al.l l2006h . 
Si/H=1.7xl0- 8 , Na/H=2.3xl0" 9 and Fe/H=1.7x 10" 9 . 

We first investigated the role of gas thermodynamics in 
the HCO + deuterium fractionation. To do this, we stopped 
to solve the thermal balance when the FUV absorption was 
large enough so that the temperature reaches a minimum 
value that we kept constant in the most shielded regions 
of the PDR. Figure [3] shows the predicted temperature 
profiles as well as the [H 2 D+]/[H+] and [DCO+]/[HCO+] 
abundance ratios as a function of the cloud depth for a 
minimum value of Tfc = 15, 20, 30 and 60 K. The predicted 
[DCO+]/[HCO+] ratio scales with the [H 2 D+]/[H+] ratio, 
as expected when DCO + gets fractionated by the reac- 
tion of CO with H 2 D + and mainly destro yed by disso- 
ciative recombination with electrons (see e.g. iGuelin et al.l 



2 [DCO+] = n(DCO + )/n(H 2 ). 
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Fig. 3. Chemical models for different minimum gas 
temperatures: 15, 20, 30 and 60 K. The density is 
?ih=4x 10 5 cm _3 and the illuminating radiation field is \ — 
60. Temperature profiles and predicted [H2D + ]/[Hj] and 
[DCO + ]/[HCO + ] abundance ratios are shown as a function 
of Ay. The [DCO + ]/[HCO + ] ratios inferred from observa- 
tions in the cold condensation at Sx ~ 40 — 45" and in the 
warm PDR gas at Sx ~ 10 — 15" are shown respectively 
with the blue and red arrows. 



1982). The displayed models also imply that low gas tem- 
peratures (< 20 K) are needed to reproduce the observed 
[DCO+]/[HCO+] ratio at Sx =40-45" (A v - 10 - 20 de- 
pending on the assumed density profile). This is easily un- 
derstood because the exchange reaction bet ween Hf and 
HP i s most efficient at low temperatures (jGerlich et al.l 
2002). Therefore, the observed [DCO+]/[HCO+] > 0.02 
abundance ratio can be reproduced using gas phase chem- 
istry only if the gas cools down from the photodissociation 
front to < 20 K, in good agreement with the DCO + ex- 
citation calculations. Note, however, that CO freeze-out is 
believed to further enhance the [DCO + ]/[HCO + ] ratio over 
the values predicted by pure gas phase fra ctionation by in- 
creasing the abundance of H t and H2D + (jBrown fc Millar! 
1989; ICaselli et al.lll999l) . The C 18 J=2-l emission shown 
in Figure [2] substantially decreases at the DCO + peak. This 
behavior is reminiscent of CO depletion but it could also 
come from a combination of lower excitation and of opac- 
ity effects. Future observations of molecular tracers of gas 
depletion are needed to constrain the dominant scenario. 

The DCO + lines stay undetected in the warm gas where 
HCO+ (not shown here) and H 13 CO+ still emit. Indeed, 
DCO + can not be abundant in the photodissociation front, 
where the large photoelectric heating rate implies warm 
temperatures (T/. > 50 K), because the reaction of H2D+ 
with H2 dominates and implies a low H2D + abundance 
([H 2 D+]/[Hj] ~ 2 x 10~ 4 ). From the upper limit of the 
DCO + emission at Sx =10-15" (Ay ~ 1), we estimate a 
low abundance ratio [DCO+]/[HCO+] < 10" 3 in the FUV 
photodominated gas, in agreement with the model predic- 
tions. 

The small distance to the Horsehead nebula (~ 400 pc), 
its low FUV illumination and its high gas density imply that 
many physical and chemical processes, with typical gradi- 
ent lengthscales ranging between 1" and 10", can be probed 
in a small field-of-view (less than 50"). The Horsehead edge 



thus offers the opportunity to study in great detail the tran- 
sition from the warmest gas, dominated by photodissocia- 
tion processes and photoelectric heating, to the coldest and 
shielded gas where strong deuterium fractionation is taking 
place. Therefore, the Horsehead edge is the kin d of source 
neede d to serve as a reference for PDR models ()Petv et al.l 
120061 ) and offers a realistic template to analyze more com- 
plex galactic or extragalactic sources. 
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